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Abstract. Wc investigate the thermodynamics of a combined Dickc- and Ising- 
model which exhibits a rich phenomenology arising from the second order and 
quantum phase transitions from the respective models. The partition function is 
calculated using mean field theory, and the free energy is analyzed in detail to 
determine the complete phase diagram for the system. The analysis reveals both 
first- and second-order Dicke phase transitions into a super-radiant state, and the 
cavity mean-field in this regime acts as an effective magnetic field, which restricts 
the Ising chain dynamics to parameter ranges away from the Ising phase transition. 
Physical systems with a first order phase transitions are natural candidates for 
metrology and calibration purposes, and we apply filter theory to show that the 
sensitivity of the physical system to temperature and external fields reaches the 
l/N Heisenberg limit. 



1. Introduction 

The understanding of the remarkable and useful properties of matter in different phases 
and of the critical behaviour near phase transitions presents ongoing challenges in 
theoretical and experimental physics [1, 2, 3, 4, 5]. Entirely new types of phase 
transition phenomena become relevant with the ability to control and engineer 
microscopic interactions and systems, e.g., in cold atom experiments with spinor Bose- 
Einstein condensates and with Fermi gases [6, 7]. Indeed, a whole branch of quantum 
computing research attempts to use candidate systems for quantum computing as 
quantum simulators in which quantum gates are operated so as to simulate suitable 
inter-particle interactions [8] and in this way implement theoretical phase transition 
models in a quantum analog computer. 

Local interactions between nearest neighbours in a spin chain lead to Ising-, 
Heisenberg-, and other interaction models with phase transitions at definite values of 
the interaction strengths and external controllable parameters, such as a bias magnetic 
field. While these interactions are reliable models of, e.g., magnetic interactions in 
solids, they can also be engineered exactly among trapped atoms or ions, with the 
added experimental possibility to control the sign and magnitude of the interactions 
with laser beams and the spin temperature by optical pumping [9]. 

Atomic and optical systems also permit the engineering of interactions between a 
large number of atoms and a quantized oscillator mode. Such systems are implemented 
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in various schemes for quantum computing, where the oscillator mode is used as 
a data bus between the atomic quantum bits. Beyond a critical coupling strength 
to the oscillator and below a critical temperature the system undergoes a phase 
transition, and the thermodynamic ground state acquires a macroscopic excitation 
of the oscillator mode. This phase-transition was first discussed by Dicke [10], and 
the Dicke phase transition has since then been studied extensively [11, 12], and was 
recently observed in experiments with cold atoms in an optical cavity [13] using 
techniques similar to those described in [14]. 

It is conceivable that one can implement both the Ising and the Dicke interaction 
Hamiltonians in many different ways using atoms and cavities, nano-mechanical 
devices or collective vibrational motion of the atoms or ions. The partition function 
of such a combined Dicke- Ising model was determined recently in the thermodynamic 
limit [15], and in this article we study the phase diagram and further properties of the 
system and possible applications. 

First order phase transitions represent discontinuous changes and hence a high 
sensitivity to variations in the values of the physical parameters of the model near 
the critical point. Sensitivity of quantum systems in fundamental metrology is a 
very active research field where much research has been devoted to identify how the 
sensitivity depends on the number of particles N. The "standard limit" 1/y/N, is 
replaced by the "Heisenberg limit" l/N within different realizations of non-interacting 
particles [16, 17, 18, 19], while more rapid decrease with particle number has been 
proposed in different models of interacting particles [20, 21, 22]. In this manuscript 
we argue for a power law decrease similar to the Heisenberg limit for measurements 
of temperature with our interacting system. 

The manuscript is organized as follows: In section 2 we outline the model and 
discuss the two limiting Ising and Dicke regimes and their phenomenology. In section 
3 we review the phase diagram calculation [15] and we recast this calculation in terms 
of mean-field theory. In section 4 we discuss the phase diagram in detail. In Section V 
we analyze the application of the phase transition for metrology. Finally, we conclude 
in section 6. 

2. The Dicke-Ising model 

2.1. The Ising Model 

Consider a one-dimensional chain of N spins or two-level atoms realizing an Ising 
model with a transverse magnetic field, 



where ft > is the transverse field and J > is the interaction strength between 
neighbouring spins and cf v n+1 = o~\. 

The model is not trivially easy to diagonalize, as the transverse field operators 
erf do not commute with the interaction terms a\o~\ +1 . In general, for h ^> J the 
transverse field dominates in which case we expect (erf) ~ and (erf) = 1. For a 
strong coupling J > /i the interaction term favours parallel spin in the y-direction 
so (erf) ?« and the system possesses two states of equal energy with (erf) = +1, or 



N 




(1) 
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An (almost) exact diagonalization of the Ising Hamiltonian can be performed 
[23, 24] using a Jordan- Wigner transformation which maps the spin-operators to 
fermionic operators. The particular mapping we employ here is c n — iY[j<n a j a n> 
where = (a* — icr^)/2 is the n'th spin-z step- up operator. Multiplying the spin- 
operators with phase factors, crj, depending on the spins at previous locations in 

the spin-chain leads to the fermionic anti-commutator relations ^c^c^/^ — 8 n , n >. 

The fermionic number operator is c\c n = CnC^ = (1 — o- z n )/2 and the state with 
all spins pointing in the ^-direction is mapped to the fermionic vacuum while S n 
flips a spin from up to down thereby creating a fermion. The inverse mapping reads 
erjj = — i rij<n(l~ 2cjc„)c„ such that the y-y-interaction term in the Ising Hamiltonian 
becomes — J erf cr? +1 = — J{c\ — c,)(cj +1 +Ci+i) for i < N. Due to the periodic boundary 
conditions we cannot remove the intermediate (1 — 2cjcj)-factors in the final y-y-term 
— Jo v N o\. This term can be written as J(cjy — cat)(c{ + c\) IljLi(l — ^ c ] c j) an d by 
adding and subtracting J(c^ N — cjy)(c\ + a) we can write the resulting Hamiltonian 
as 

Rising = - 2c\ Ci ) -J^2 C l Ct +! + C Ul C i + C »+1 C » + 44+1 

i i 
N 

+J(ci - c N )(c\ + Cl )(l + 2c] Cj )) 

where the last term is of relative order 1/N since 1 + Yij(^ ~ 2cjcj) is a projection 
operator onto the subspace of even total spin in the ^-direction and therefore is of 
order 1. 

By neglecting this less important final term the Hamiltonian is quadratic in 
the fermionic operators, and can be diagonalized using a Bogoliubov-transformation. 
The Bogoliubov-transformation can be decomposed into a single-particle basis change 
to momentum space and quasi-particle operators connecting particles of opposite 
momenta [23, 24, 25] and we get: 

Rising = £ e(k) U lk 0(l/N), (2) 

where the operators jk are the Bogoliubov transformed fermionic operators, e(fc) = 
2( J 2 + h 2 — 2Jh cos(fc)) 1 / 2 and k e has the form 2im/N and is a reciprocal lattice- 
vector to the chain-lattice. The eigenmodes of -fusing correspond to free fermions with 
a dispersion relation given by e(k). Note that (2) has a single unique ground state in 
contrast to (1) which has a two- fold degeneracy. This difference is due to the omitted 
term and will not be relevant in the thermodynamic limit. 

Since the spectrum of the Ising model is so simple we can calculate the partition 
function for the system. The partition function is given by Z® sing = Tr(cxp(— /3i? Ising )) 
where f3 is the inverse temperature. We can easily calculate this trace, 

Z° ing (/3, h, J) = [] Tr [cM-Pmilllk - 1/2))" 

fceZjv 

= Y[ 2cosh( / 3e(fc)/2). 

fceZjv 

Using this result we can also calculate the free energy given by i*i s ; ng = 
— /3" 1 logZj° sing . The diagonalization of i?i s i ng is exact to order 1/N and to the same 
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precision we can replace the sum over k e Z^v by an integral, 



F Uins (l3,h,J) = ~ log(2cosh(^e(fc)/2)) 



~-N —log(2cosh(pe(k,h,J))). 



(3) 



J —71 



Note that /3e(fc, h, J)/2 = (3J(l + (h/J) 2 -2(h/J) cos(fc)) 1 / 2 , and hence it is convenient 
to parametrize Fi s i ng with the dimensionless quantities (3 — (3J and h = h/J in the 
following. 

The Ising model exhibits an infinite order quantum phase-transition between a 
paramagnetic h^> J and ferromagnetic J ^> h phase with critical point at J = h where 
a non-analyticity arises in e(k) at k = 0. This system has been studied extensively [25] 
and is a prime example of a quantum phase transition. The quantum phase-transition 
does not survive to finite temperatures, but it has important consequences for the 
finite-temperature behaviour of the system near the critical point. 

2.2. The Dicke model 

Another phase transition model consists of spins or two-level atoms coupled to a 
harmonic oscillator mode with frequency w, 



where of is the Pauli x-matrix acting on the i'th spin and a is the harmonic oscillator 
step down operator. We denote the coupling-strength of the oscillator to a single spin 



scaling thus corresponds to atoms with a constant spatial density which is well defined, 
also in the thermodynamic limit N — > oo. 

The explicit form of the interaction is well established in quantum optical systems, 
and occurs both for two-level systems and for, e.g., Raman processes between two 
states via an intermediate excited state through absorption and stimulated emission 
of the quantum field and a classical control field. We note that in the Jaynes- 
Cummings model of a single two-level atom and a single field mode, the rotating wave 
approximation retains only the terms a^a + oa) in V and provides a considerable 
simplification of the problem. The partition function for the Dickc-modcl has been 
calculated analytically in the thermodynamic limit and a second order phase transition 
has been identified [11]. While that calculation pertained to the rotating wave 
approximation it has been shown [12] that this does not change the Dicke phase- 
transition qualitatively. Within our application of a mean field approximation to the 
combined Dickc-Ising model we shall retain the full interaction (5) as, the rotating 
wave approximation is more difficult to deal with. 

The Dicke model H osc + V can be realized experimentally as described in [14] 
where a dynamical version of the standard Dicke model is investigated in a cavity 
using four-level atoms coupled by Raman channels. The model parameters of this 
system are functions of the atomic and field parameters applied and can be tuned 
over large ranges. 



H osc = u>a'a, 
through the interaction 



(4) 




(5) 




Phase transitions in an Ising chain interacting with a single mode cavity field 



5 



3. The partition function for the Dicke-Ising model 

The combined Dicke-Ising model has total Hamiltonian H = -ffising + -Hose + V with 
three parameters g, h, and J which can be varied independently in atomic simulators 
of the model. With J = the Hamiltonian realizes the Dicke model of N two- 
level atoms interacting with a cavity-field via the dipole interaction. In this regime 
2h correspond to the energy-splitting of the individual two-level atoms and g is the 
coupling strength to the cavity-field. 



3.1. Coherent state integral 



We will proceed by writing first the expression of the partition function for the full 
problem Z = Tr(cxp(-/3 J ff)) with H = i?i sing + H osc + V. Following [11] wc write H 

as 



JV 



-h*f-JoVoV +1 +g*l 



it 



N 



it 



N VN 



We will evaluate the trace over the oscillator-mode in the partition function using 
the coherent state representation of the field, and for this purpose we should bring 
exp(— j3H) on normal ordered form with respect to the scaled operators b = a/y/N 
and = a) /\fN. Due to the scaling, their commutator is [b,tf] = 1/N and vanishes 
for N — > oo. 

The approximation of neglecting the commutator in the thermodynamic limit 
can be illuminated by Wick's theorem [26]: Any string of creation- and annihilation 
operators can be written as their normally ordered form plus additional terms. These 
terms are given by means of a contraction defined by C(AB) = AB—-.AB:, where 
: AB: is the normal ordering of AB. Using this notation, Wick's theorem states that 

ABC ...Q =:ABC ...Q: + J2 : 

one contraction: 



two contractions: 



Now, in our case C(btf) = btf - tfb = [b,tf] = 1/N and C(tfb) = 0. Hence, 
applying Wick's theorem to the expansion of the exponential in Z = Tr(exp(— f3H)), 
and assuming that the limits in Z = limjv->oo li m i?^oo J2^=o(~ PH) r /rl can be 
interchanged, we see that all terms involving contractions will be of order 1/N or 
higher and can therefore be neglected in the thermodynamic limit. More generally any 
expression of the form Tr(/(6, tf)e~@ H ), where / is a polynomial, can be calculated 
to an accuracy of 1/N using the normal order : f(b, b^)e~^ H :. The validity of this 
truncation for the Dicke-model was discussed in [12]. 

When performing the trace of a normally ordered operator it is convenient to use 
the coherent states, \a), eigenstates of the annihilation operator: a\a) = a\a), which 
yields 



-/3u|ap 



Tr, 



spin 



exp 



-haf — J a 



i+i 



+ 



2g^{a) 



N 



The remaining trace over the spin degrees of freedom is exactly equivalent to the the 
original Ising model calculation, where the term involving the real part of the complex 
field argument acts as an additional magnetic field in the x direction, and hence the 
Ising model is biased by an effective magnetic field in the xz-plane with magnitude 



+ 0(1/N). 
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h 2 s = h 2 + 4g 2 3ff(a) 2 /A r . Since this effective field lies within the plane orthogonal to 
the t/-axis we can choose the direction of the effective field as a redefined z-axis and 
apply the same diagonalization as in the pure Ising problem. 

This yields Z = J ^ exp(— f3u \a\ 2 + log Zj° sing (/3, h c g, J)) and after performing 
the substitution z — a/y/N we get 

Z = N J ^ exp (-Npu |z| 2 ) 

•exp ^ y" r dfclog(2cosh(/3e(fc,/ lcff (5R(z)), J)/2))j . (6) 

where the dispersion relation e(k, h e s(x), J) is for the effective magnetic field h e s and 
spin coupling J. 



3.2. Mean field theory 

We may also attack the original problem with an Ansatz replacing the interaction part 
V of the Hamiltonian by the mean field expression 



Vmf = gNs x (a + a!)/VN + 2g ^ afx - 2gNs x x 

i 

where the c-number mean fields read s x — (J2iLi a f) /N and x = (a + a)) /2y/N. The 
last term compensates for double-counting of the interaction energy, while correlations 
in the mean-field fluctuations g(^2 { erf — Ns x )(a + — \/Nx)/\/N are omitted. 

With the Hamiltonian in this form the mean-fields split the Hamiltonian into two 
separate terms: a classically driven field mode -/Jmf.osc = wa'a + gNs x (a + a))/\fN 
and an Ising model with an additional magnetic field component along the x-direction, 
#MF,isin g = —h °~i + %9 X ~ J J2i a i a i+i- Tnis mean field Ising Hamiltonian 

has a field in the xz-pl&ne with magnitude ft. 2 ff = h 2 + 4g 2 x 2 . The partition function 
for the mean field Hamiltonian therefore factors Z = Zmf.osc-^mf, ising exp(2(3gNs x x) 
and can be readily determined for arbitrary values of the mean field amplitudes x and 
s x . The corresponding free energy reads 

F = \ log(l - e-^) - - 2Ngs x x + Nf lsing ((3, h eS , J), 

p CJ 

where /i s i ng = Fi s i ng /N is the free energy per particle for the atoms as calculated in Eq. 
(3). Note that the free energy includes a term representing the thermal distribution of 
the cavity-photons as well as the contribution depending on the field amplitude. The 
values of the mean-fields can now be obtained by minimizing the free energy. A short 
calculation reveals 

1 dF 2g 2 s x 



N ds T co 



— 2gx (7a) 



1 9F dj ' I sing 

Ndx-=- 2gS * + ^x~> (7b) 
This shows that s x = —xiojg and that the Dicke-order parameter x should be found 
by minimizing cox 2 + /ising with respect to x. We will return to this minimization 
problem below. 

In the mean-field description the interpretation of the physical properties of the 
system becomes clear and unambiguous. As an example, with the mean-field theory 
we can obtain expressions for various correlation functions for the atomic variables 
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x 2 for J + g + h = 1.0, UJ = 1.00, = 100.0 % for J + g + h = 1.0, U = 1.00, /3 = 100.0 




Figure 1. (Colour online) Observables for the Dicke-Ising model with = 100, 
lj = 1 and g + J + h = 1. Left: order parameter x 2 and right: susceptibility along 
the z-axis \ = dm/dh. All dotted lines converging to a corner correspond to a 
fixed ratio of the parameters indicated on the edge opposing the corner. The black 
solid line is where the order parameter x becomes non-zero, corresponding to the 
Dicke phase transition for J = 0, and the red solid line is where the h a ff/J = 1, 
corresponding for g = to the Ising phase transition. 



from the large amount of theory already present on the transverse Ising chain since 
in thermodynamic equilibrium the strongly coupled system is effectively identical to 
a rotated Ising chain. 

It is reassuring, but hardly surprising, that the mean field result can be recovered 
from the coherent state integral for the thermodynamic limit N — > oo result. By 
Laplace's (saddle point) method, one can replace the integral over coherent state 
amplitudes by discrete contributions from the location of the maximum of the 
exponential with respect to z in (6). It is easily shown that this maximization coincides 
with the mean field, identified by minimization of the free energy and by picking one 
of the maximizers in (6) to represent a symmetry-broken physical state of the system. 

4. Phase diagram and analysis of the free energy 

We now turn to the problem of finding the minimum of the free energy (7a), (76) along 
with several important observables like the magnitude of the oscillator mean field, the 
spin magnetization and the susceptibility \ — dm z /dh. 

For fixed j3 and w, the system is controlled by three parameters h, J and g. To 
illustrate the phase transitions in the system, we introduce a convenient way to plot 
different quantities as function of these variables in Fig. 1 and 2. In each plot, the 
sum h + J + g is fixed, and the corners of the triangles shown correspond to each 
of the three quantities acquiring the maximum value while the others vanish. The 
straight dotted lines converging to the corners of the triangles correspond to definite 
values of the ratio between the two quantities indicated on the edges of the triangles. 
These plots can be thought of as slices of the three-dimensional simplex defined by 
h + J + g + uj = e/(3 where e should be interpreted as the system energy-scale. In this 
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Figure 2. (Colour online) Observables for the Dicke-Ising model with ,3 = 4, 
uj = 0.62 and g + J + h = 0.6. Left: order parameter x 2 and right: susceptibility 
along the z-axis \ = dm/dh. All dotted lines converging to a corner correspond 
to a fixed ratio of the parameters indicated on the edge opposing the corner. The 
black solid line is where the order parameter x becomes non-zero, corresponding to 
the Dicke phase transition for J = 0, and the red solid line is where the h c g/ J = 1, 
corresponding for g = to the Ising phase transition. 



coordinate representation, we show with colour coding the value of different interesting 
quantities. 

Of particular interest in Fig. 1 is the oscillator field strength, represented by x 2 , 
and the susceptibility \ which are shown for the case where w = 1, /3 = 100 and 
h + J + g = 1 for Fig. 1 and w = 0.25, /3 = 4.0 and h + J + g = 1 for Fig. 2 where h, 
J and g are positive. 

The edge of the triangle between g and h (i.e., with vanishing J) corresponds 
to the usual Dicke- model, while the edge between h and J (i.e., with vanishing g) 
corresponds to the usual Ising-model with the critical point at h/ J = 1, showing up 
clearly as a signature in the variation of the susceptibility \- We observe that this 
signature is present also for finite Dicke coupling parameter in the plot. The black 
curve in each plot shows where the Dicke phase transition occurs. In both Fig. 1 
and Fig. 2 signatures of both first- and second-order phase transitions can be seen. 
Approximately below the line h/J = l the second order transition can be identified 
by the smooth increase in x 2 whereas above h/J=l one can discern a discontinuous 
jump in the order parameter. 

The most significant difference between the cases presented in Figs. 1 and 2 is 
the susceptibility \. For moderately low temperatures, j3 = 4 in Fig. 2, the signature 
of the Ising quantum phase transition is still clearly present, whereas for very low 
temperatures, /3 — 100 Fig. 1, the Ising phase transition becomes almost completely 
suppressed in the Dicke regime. Indeed, far into the Dicke regime (towards the right 
vertex in the triangles) the spin interactions do not appear to play any significant role. 

Looking at the black and red solid lines in Fig. 1 one might be tempted to conclude 
that in the super-radiant phase we always have h e g/J > 1 (i.e. the red line does not 
penetrate into the area to the right of the black line). A close look at Fig. 2 will, 
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however, reveal that this is not always the case. Indeed, for h + J + g = 0.6, uj = 0.62 
and (3 = 4 there is a small part of the parameter-space where x > and h c g < J. 
The physical reason why the Dicke phase transition almost, but not quite, suppresses 
the Ising transition remains to be understood. Formally, it occurs because the free 
energy has a minimum giving an effective magnetic field such that h e g/J > 1. The 
Ising critical point is therefore simply skipped in these cases and only an amputated 
signature of the Ising phase transition is present in cases such as shown in Fig. 1. 

To investigate this phase transition in more detail, let us further consider the 
minimization of the free energy. The order-parameter x enters the mean field Ising- 
term via the effective magnetic field (h e s/J) 2 = h 2 + 4g 2 x 2 /J 2 . If we therefore 
introduce a rescaled order-parameter x — 2gx/J and a rescaled mode frequency 
oj = ujJ/Ag 2 we can write the free energy as a function of a few dimensionless quantities 

^TF = - T -^log(2cosh(^(l + (h 2 + x 2 ) - 2(h 2 + x 2 ) 1 ' 2 cos(k)) 1 / 

+LOX 2 + C, 

where C is a constant independent of x. 
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Figure 3. (Colour online) Examples of various functional shapes of (_F — 
F(0))/NJ as a function of h e g. The different line styles indicate varying Cj and 
the two colours indicates P = 1 and P = 10 respectively. 

In order to understand the structure of the phase diagram it is necessary to 
investigate how the integral changes as a function of j3 = f3/J and h compared to uix 2 . 
If we choose the variable h a s = h e g/ J = \/h 2 + 4g 2 x 2 / J as the independent variable 
instead of i, the integrand in (8) only depends on (3 and the new variable h c g, while 
h e ff > h imposes a boundary condition on the minimization with respect to h e g. To 
avoid confusion, we will consider F a function of i and use the symbol F to denote the 
dependence on h e g. The system is in the super-radiant phase whenever the minimum 
in F occurs for h c s > h which implies i^0. Examples of F for representative values 
of Co and /? can be seen in figure 3. By a numerical investigation it is quickly revealed, 
that F has at most a single local minimum (e.g. the curves (a), (b), (c) and (d) in 
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figure 3) at h e ff ^ or no local minimum (curve (e) in figure 3). The existence and 
location of the minimum are thus solely determined by /3 and Co. 

This implies that when keeping oj and (3 fixed the minimum of the free energy is 
either at h c s = h or at the local minimum of F. If we imagine tuning h from high 
values towards low values (i.e, setting the boundary condition h c g > h at different 
locations, for example along curve (c) in figure 3) the system will pass a second order 
phase-transition when h passes the local minimum of F. In the case that the local 
minimum is not the global minimum there will be a h > where F goes below its value 
at the local minimum implying that when one further lowers h the system will undergo 
a first-order phase transition into the normal state again. By the same reasoning if F 
has a single global minimum (curve (a) and (b) in figure 3) there can only be a second 
order phase transition when tuning h. 

The second-order phase transitions can be investigated in further detail using 
Gintzburg-Landau theory: In the neighbourhood of the second-order phase transition 
the order parameter x will always be small so we can expand the free energy F as a 
polynomial in x around x = 0: 

^1^C + I o 0, h) + (Q + I 2 (fi, h))~x 2 + I 4 0, h)x\ (9) 

where I n is the n'th term in the Taylor expansion of the integral (8) with respect to x. 
The standard argument from Gintzburg-Landau theory is now that this fourth-order 
polynomial has a non-zero minimum when uj + I 2 (f3, h) is negative. The second order 
phase-transition therefore occurs when I 2 {j3, h) = — <D. By numerical investigation one 
finds that — I 2 is bounded by approximately 0.3356 implying that for u> > 0.3356 no 
phase transition can occur. 

The first-order phase transition grows out of the second-order phase transition so 
there will be a region where the first-order jump in the order parameter is small. In 
that case we can still use Ginzburg-Landau theory and in particular we can find the 
point where the second-order transition changes to a first-order transition, i.e. when 
a local minimum in F changes from purely local to truly global. Again we analyze 
the polynomial expansion and one can show [27] that one needs to solve the system 
of equations I ±{13, h) = and h) + Co = to obtain the point where the phase 
transition changes nature. By investigating the functional form of I4 it turns out that 
there is a minimal f3 c below which the first order phase-transition cannot occur. This 
value can be calculated numerically and is approximately /3 C w 1.1430. 

This identifies where the second-order transition changes to a first-order 
transition. To determine the first-order transition boundary for finite jumps in the 
order parameter, however, it is necessary to deal with the free energy F to all orders. 
Numerically it is not difficult to investigate for which value of h the value of F coincides 
with the value at the local minimum as described above. All this information has been 
combined into figure 4 where the phase-boundaries for various values of u> have been 
indicated. The coloured dotted lines represent first order transitions whereas the 
solid lines indicate second-order transitions. The black dashed curve indicates where 
Ii{(3,h) = and its intercept with the curves I 2 {(3,h) + u> = indicates where the 
phase transition changes type between first- and second-order transitions. 

5. Using a phase transition for a high precision measurement 

A first order phase transition is interesting for many different reasons and here we 
consider its use as a measurement tool. Indeed, the standard description of a first 
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Figure 4. (Colour online) Phase diagram for the Dicke-Ising model. Each curve 
shows the location of the phase transition for the specified parameters. If the 
curve is full the transition is second order, whereas a dotted curve represents a 
first order transition. The dashed black curve represents the general boundary 
h) = where the transition changes from second to first order along each of 
the coloured curves. 

order phase transition includes a discontinuous jump in the order parameter, and it 
is a relevant question, how precisely an experiment can locate the position of this 
discontinuous jump. The size of the jump-discontinuity usually scales linearly with 
the number of particles, while the width of the transition region often scales with an 
inverse power of this number, and under that assumption we shall present a simple 
model for the metrological sensitivity of the system. Since the phase transition occurs 
for rather non-trivial combinations of the temperature and the interaction parameters 
of the models, by changing some of these parameters in a controllable way, one may 
be able to select parameter ranges with particularly high sensitivity of the phase 
transition point to the value of the quantity being probed. 

We will consider a measurement strategy where the system is probed at some 
range of values of some control parameter, e.g. a bias magnetic field. Fig. 5 shows 
how the cavity field order parameter varies as a function of h for three different 
temperatures. To produce this plot, we have selected values of Cj and thus an 
area of the phase diagram where the critical magnetic field depends strongly on the 
temperature, cf. the steepness of the dashed curves in Fig. 4. We expect that it will be 
possible to determine the critical value of the magnetic field with high precision, and 
since in this case a variation of the temperature of 1% changes the critical magnetic 
field by approximately 35%, the measurement of the critical field yields a very sensitive 
temperature measurement within the appropriate range of values f3 « 0.77 ± 1%. 
Sensitivity in, e.g., a lower temperature range is obtained if we chose a higher value 
of w and scan a different range of values of the magnetic field. 

The order parameter presented here is the intra cavity field intensity. We imagine 
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Figure 5. (Colour online) The order-parameter function of h and for 

the lower-right part of the curve for lj = 0.301 in Fig. 4 where the slope of the 
phase transition line is large. For the blue curve 1/0 = 0.77, for the green curve 
the temperature is 1% higher and for the red curve the temperature is 1% lower. 

that the cavity leaks photons at a sufficiently low rate not to significantly disturb the 
thermodynamic steady state of the system, and herewith, detection of the intensity 
of the emitted light is a direct probe of the cavity field order parameter. We assume 
that the inverse temperature /3 is known to be close to some reference value, and we 
can then estimate the difference 5/3 by the best unbiased linear estimator as described 
in Appendix Appendix A. This estimator is given by 

where rij is the detected number of photons in a given time while the controllable 
effective bias field h attains the value hi- fJ,(hi) is the expectation value of the photon 
number and o~ 2 (hi) is the photon number variance. /jf(h) denotes the derivative of 
the expected photon number with respect to changes in inverse temperature /3, and 
the expression applies within a narrow range where a linear variation of the expected 
photon number with (3 is valid. 

In the limit of high bias field resolution, the sum in the estimator can be 
converted into an integral, and one can determine the variance of the estimate: 

V&r(5(3(n(h)) = 1/ J fj,' (h) 2 /cr 2 (h)ih, see details in the appendix. 

So far the arguments have been of a general nature. Let us now assume the Dicke- 
ring model, in which the photon number distribution is well described as a thermal 
state below and a displaced thermal state above the Dicke phase transition. The first 
and second moments of such distributions can be calculated using, e.g. the positive 
P-representation for the thermal state, 

H(jk,p) = (n) = Nx 2 + n 
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o- 2 (h,f3) = Vara, = Nx 2 (l + 2n) + fi + n 2 

Recall that the order parameter x 2 is a function of the system parameters h, P and 
Q and in the thermodynamic limit it has a discontinuous jump at the dashed lines 
shown in Fig. 4. For a finite system, however, the phase transition constitutes a 
smooth curve with a fast increase of the order parameter. The width of this region 
is not easy to determine but finite size effects in phase transitions tend to smoothen 
phase transitions leading to a decreasing width as N increases. Indeed, [28] finds that, 
in general, the width scales as 1/N, but for the sake of generality we assume a scaling 
JV-t, 7 > 0. 

Since both (n) and Vara are proportional to Nx 2 both n(h), ^'(h) and a 2 (h) 
will carry a signature of this power law. To be explicit assume that Nx 2 {h,P) = 
7Vtanh(iV 7 (/i — h c (/3)) where h c is the critical value of h c as a function of /?. Then 
fj,(h) cx N, n'{h) = (dp(j,)(h, fa) cx N 1+ ^h' c o ) and a 2 {h) cx N. The function // only 
has support in a region of width 1/N 1 near ho = h c {Po)- The variance of the estimate 
6p then scales as V&rdp w ((/i / (ft ) 2 /<7 2 (ft ))JV-T)~ 1 cx (h' c (M 2 N 2+2 '< / N)- 1 N 1 , i.e. 

1 

Var<5/3 cx - 



h' c (p ) 2 N^ 

This is our main result of this section, showing that the sensitivity is better than 
the "standard limit" where one expects Var/3 <~ 1/N, and depending on the character 
of the finite size effects (the power 7) , it is potentially also better than the Heisenberg 
detection limit. With the result from [28], 7=1 the accuracy is actually at the 
Heisenberg limit. Note that the above argument is quite general and applies to any 
first order transition with an intensive order parameter. 

The term h' c is included in our expression in order to show explicitly that the 
sensitivity depends on the curve of critical points in the phase diagram of the system. 
From Fig. 4 we see that h' c can be chosen large for arbitrarily small temperatures by 
tuning Co. The large value of h' c , however, comes at a cost: The slope of h c is highest 
near the thick dashed curve, which is also where the first order transition has small 
amplitude and changes to a second order transition. In a concrete implementation, 
the values of Cj and the range of effective magnetic fields need to be chosen with care 
to reflect the actual scaling N~ J and the size of the jump discontinuity. 

With an adaptive measurement scheme, we imagine that the number of iterations 
with different hi for a reliable detection of the critical value of the /i-parameter can 
be optimized. It is clear that a more detailed investigation is necessary in order to 
quantify the accuracy and scaling of resources of such measurements. Indeed, the 
specific power law 1/N 1 for the transition width is only a convenient Ansatz, and a 
non-mean field calculation on a finite system will be needed in order to investigate the 
approach towards the thermodynamic limit in more detail. Furthermore, the critical 
properties and the long range correlations of the system may possibly lead to even 
better estimates by use of the recent techniques of quantum non-linear parameter 
estimation [29, 21]. These issues we shall defer to a later publication. 



6. Conclusion 



We have investigated the thermodynamic properties of a Dicke-Ising model 
incorporating both the quantum transverse Ising model and the Dicke model as special 
limiting cases. We have derived expressions for the free energy using a coherent state 
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integral similar to [11] but also using a mean field theory with a clearer interpretation 
for the field statistics. The combined model exhibits a first order phase transition 
which is not present in either of the two separate models. By a simple numerical 
search the free energy minimum can be identified and the value of the Dicke mean 
field and the magnetic susceptibility can be determined as functions of all physical 
parameters of the model, cf. Figs. 1 and 2. 

Using the free energy and Ginzburg-Landau theory we also investigated the 
complete phase diagram as shown in Fig. 4. The Dicke phase transition occurs 
also for moderate inter-particle interactions and the Ising phase transition is also 
well preserved for weak and moderate light-matter couplings. In a small area of the 
parameter space both phase transitions coexist closely together, but for a stronger 
Dicke model interaction the resulting mean field puts the system in a regime without 
any observable Ising phase transition. 

The Dicke- Ising model constitutes an interesting mix of second-, first- and infinite 
order phase transitions. The interplay of these phase transitions and a complete 
description beyond the mean field approximation of the fundamental excitations at 
the critical points would be and interesting continuation along the lines of this work. 
In addition to its fundamental theoretical interest, the first order phase transitions 
provides a tool for precise measurements of, e.g. the magnetic bias field or of 
the temperature. We have presented a simple estimate of the accuracy of such a 
measurement device showing that the variance scales as 1/N 1+1 which is better than 
the standard limit 1/N for independent measurements on TV particles. 

This work was supported by the European Union Integrated Project AQUTE. 

Appendix A. Best unbiased linear estimator 

We will consider the best unbiased linear estimator in a situation where an 
experimenter performs a sequence of measurements where she scans a parameter x 
(e.g. a bias magnetic field) in order to uncover another, unknown, parameter q (e.g. 
the location of a critical point). The experimenter measures a discrete stochastic 
variable n (e.g. photon number) which has probability distribution p(n;x,q). In the 
measurement the experimenter thus collects, for a fixed q, the values n, corresponding 
to selected Xi . 

In the following we assume that the unknown parameter q is close to a reference 
value, which, without loss of generality, we take to be zero. For g C 1 we can then 
expand the moments of p in a Taylor expansion such that 



where we have expanded the mean p,{x,q) to first order (and p,(x) = fi(x,q), 
fi'{x) = (d q n)(x, 0)) and the variance cr 2 (x, q) to zeroth order in q. 
A linear estimator is of the form 



where {<?i} are weighting coefficients, {n^} are the observed values of rij corresponding 
to the chosen values Xi and c is a constant. To find an unbiased estimator we require 
E [q({ni})] cx q which implies c = — J2i To find the best linear estimator we 

optimize the signal-to-noise ratio E [<?({^i})] 2 /Var(g({ni})) with respect to the vector 



E[rii] = fj,(x,q) « fi{x) + qn'{x) 
Var(rii) = a 2 (x 7 q) rts o- 2 (x) 
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Qi. To second order in q the signal-to-noise ratio is q 2 (^2 i gifjl (%i)) 2 / X)i 9i a2 { x i)- The 
minimum of the signal-to- noise ratio is obtained when gi oc /i' (xi)/a 2 (xi). 

The constant of proportionality should then be chosen such that q = E [q({ni})] = 
^2 i gi(E [rii] — n{xi)) which gives the condition ^\ 9iH'( x i) — 1 an d the normalization 
constant is given by A^ 1 — ^ii^'i) 2 / a f ■ With this normalization the variance of the 
estimate q is 

Var^K})) = ^ { J^pj (A.l) 
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